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Light Quark Masses in QCD 

Eduardo de RafaeP * 

'^Centre de Physique Theorique, CNRS-Luminy, Case 907, F-13288 Marseille Cedex 9, France 

I report on recent work done in collaboration with Laurent Lellouch and Josep Taron |lj concerning the 
derivation of rigorous lower bounds for the combination of light quark masses tUs + m„ and nid + mu- 
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1. INTRODUCTION 

The light u, d and s quarks in the QCD La- 
grangian have bilinear couplings 



*-'QCD — '-'QCD 



^ mgq{x)q{x) , 

q—u,d.s 



(1) 



with masses ruq ^ which explicitly break the 
chiral symmetry of the Lagrangian. Twenty five 
years of efforts of many people (including signifi- 
cant contributions from the chairman of this Con- 
ference Stephan Narison,) have resulted in the fol- 
lowing values 1^: 

TO„(4GeV^) + mrf(4GeV^) = (9.5 ± 1.9)MeV, (2) 
and §1 

TO„(4GeV2) + m,(4GcV2) = (149 ± 25)MeV. (3) 

These are values of the running masses at a scale 
ji^ = 4GeV^ in the MS renormalization scheme. 
[The choice of scale is of course arbitrary; 4GeV^ 
has become the conventional choice of the lattice 
community. We shall adopt it here as well so 
as to make easier the comparison of our bounds 
with their results.] There are quite a few inde- 
pendent determinations, including some from nu- 
merical simulations of lattice QCD , which now 
agree within errors with these values. 

The values in (||) and @) are smaller than those 
obtained in the very early determinations of the 
light quark masses using QCD sum rules (see e.g. 
ref. H] and references therein.) This is partly due 
to the fact that the experimental determination 
of ■'^jjjg has resulted in larger values than those 

*Invited talk given at the conference "QCD 97", Montpel- 
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quoted in the early 80's, and partly to the ap- 
pearance of rather large coefficients in the pertur- 
bative expansion of the relevant QCD two-point 
functions. 

Since the summer'96 there have appeared some 
determinations of the light quark masses from 
new analyses of numerical simulations of lattice 
QCD 1^,^ which find substantially lower values 
than the ones quoted above. More recently, an 
independent lattice QCD determination using 
dynamical Wilson fermions finds results which, 
for TTiu + nid are in agreement within errors with 
those of refs. 0,^; while for nis they agree rather 
well, again within errors, with the value quoted 
in (H). The present situation concerning the val- 
ues of the light quark masses has thus become 
rather confusing again, and this is why we have 
decided to concentrate on what can be said on rig- 
orous grounds about the light quark mass values 
at present, rather than argue in favour or against 
the reliability of a particular determination. 

We shall show that there exist rigorous lower 
bounds on how small the light quark masses mj -|- 
m„ and nid + can be. The derivation of the 
bounds is based on first principles: the fact that 
two-point functions 

*(q2) = i J d^xe'^-^{0\T {O{x)O{0)^) |0) (4) 

of local operators 0{x) obey dispersion relations, 
and the positivity of the corresponding hadronic 
spectral function 



x(o|o(o)|r)(r|ot(o)|o). 



(5) 
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where the sum over F is extended to all possible 
on-shell hadronic states with the quantum num- 
bers of the operator O. 

We shall be concerned with two types of two- 
point functions, one where O is the divergence of 
the strangeness changing axial current 

0{x) = df_,A'^{x) = {ms + niu) ■.s{x)ij5u{x): ; (6) 

and another one where O is the scalar operator 
S{x), defined as the isosinglet component of the 
mass term in the QCD Lagrangian 

0{x) = S{x) = m[:u{x)u{x): + ■.d{x)d{x):] , (7) 

with 

m=^{mu+md). (8) 

We shall refer to the corresponding two-point 
functions respectively as the "pseudoscalar chan- 
nel" and the "scalar channel". Both two-point 
functions obey dispersion relations which in QCD 
require two subtractions, and it is therefore ap- 
propriate to consider their second derivatives 



dt— 

(t 



^ r-Im*(t). (9) 
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The bounds follow from the restriction of the sum 
over all possible hadronic states which can contri- 
bute to the spectral function to the state (s) with 
the lowest invariant mass. It turns out that, for 
the two operators in (H) and (^, these hadronic 
contributions are well known phcnomenologically; 
either from experiment or from chiral perturba- 
tion theory (xPT) calculations. On the QCD side 
of the dispersion relation, the two-point functions 
in question where the quark masses appear as an 
overall factor, are known in the deep euclidean re- 
gion where ^ ^qcd from perturbative QCD 
(pQCD) at the four loop level. The leading non- 
perturbative ^ -power corrections which appear 
in the operator product expansion when the time 
ordered product in is evaluated in the physical 
vacuum are also known. As we shall see, the 
bounds we find cast serious doubts on the reliabil- 
ity of the recent lattice determinations reported 
in refs. 



2. PSEUDOSCALAR CHANNEL 

We shall call ^'5((?^) the two-point function in 
associated to the divergence of the strange- 
ness changing axial current in (|^). The lowest 
hadronic state which contributes to the spectral 
function in (||) is the i^-pole. From eq. (||) we 
then have 



(Mi + Q2)3- 

2 1 
dt- ^^-lm^5{t) , 



(10) 



/to 't + Q^rn 

where to = {Mk + 2771,^)^ is the threshold contin- 
uum. 

It is convenient to introduce the moments 
SAr(Q^) of the hadronic continumn integral 

2 



dt- 



t + Q^j (11) 

One is then confronted with a typical moment 
problem (see e.g. ref. The positivity of 

the continuum spectral function ^lm'i>c^{t) con- 
strains the moments Y,]y{Q'^) and hence the l.h.s. 
of ( p^ ) where the light quark masses appear. The 
most general constraints among the first three 
moments for = 0, 1, 2 are: 

Eo(Q')>0, Si(Q')>0, S2(0')>0; (12) 

Eo(Q')-Si(g2) > 0, ,Si(Q2)_S2(Q2) > 0;(13) 



So(Q')S2(Q2)-(Si(Q2))'>0. 



(14) 



The inequalities in eq. ( [i3| ) are in fact trivial un- 
less 2Q^ < to, which constraints the region in 
to too small values for pQCD to be apphcable. 
The other inequalities lead however to interest- 
ing bounds which we next discuss. 

2.1. Bounds from Eo(Q^) > 0. 

The inequality EqIQ^) > results in a first bound 

on the running masses: 

r .^9m2 167r2 2/2Af4 
[ms{Q^) + m^iQ^)] > — ^ ^4^x 



1 



1 



11 a.(Q^) 

3 TT 



(15) 
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Notice that this bound is non-trivial in the large- 
Nc hmit (/I- 0{Nc)) and in the chiral hmit 
{rUs ~ M]^). The bound is of course a function 
of the choice of the eucUdean Q-value at which 
the r.h.s. in eq. ( [l5| ) is evaluated. For the bound 
to be meaningful, the choice of (5[GeV] has to be 
made sufficiently large so that a pQCD evaluation 
of the QCD factor 



(16) 



associated to the two-point function is ap- 
plicable. The dots in (16) represent higher or- 
der terms which have been calculated up to 
0{al) [Q, as well as non-perturbative power 

corrections of I^jIiII; ^ind strange quark 

mass corrections of O y^t) a-nd O {j^^ includ- 
ing 0{as) terms [^6|-^,|j. We consider that, 
Q > 1.4 GeV is a safe choice for a numerical es- 
timate of the bounds: the evaluation of at 
Q = lA GeV with or without inclusion of the 
0{a^) terms differs by less than 10%. 

The lower bound which follows from eq. ([l5| ) for 
rUu + rUs at a renormalization scale — 4 GeV^ 
results in the solid curves shown in Fig. 1 be- 
low. These are the lower bounds obtained by let- 



all the numerical analyses of this work. Notice 
that the hadronic continuum integral in (|ll|) is 
always larger than the contribution from three 
light flavours; and indeed, we have checked that 
using Uf = 4 in the QCD expressions leads to 
even higher bounds. Values of the quark masses 
below the solid curves in Fig. 1 are forbidden by 
the bounds. For comparison, the horizontal lines 
in the figure correspond to the central values ob- 
tained by the authors of ref . ||^ : their "quenched" 
result is the horizontal upper line; their "un- 
quenched" result the horizontal lower line. 

We wish to emphasize that in Fig. 1 we are 
comparing what is meant to be a "calculation" 
of the quark masses -the horizontal lines which 
are the lattice results of ref. Q- with a bound 
which in fact can only be saturated in the limit 
where Eo(Q^) = 0. Physically, this limit corre- 
sponds to the extreme case where the hadronic 
spectral function from the continuum of states is 
totally neglected! What the plots show is that, 
even in that extreme limit, and for values of Q in 
the range 1.4 GeV < Q < 1.8 GeV the new lattice 
results are already in difhculty with the bounds. 

2.2. Bounds from the Quadratic Inequality 

The quadratic inequality in ( p^ results in im- 
proved lower bounds for the quark masses which 
we show in Fig. 2 below. Like in Fig. 1, these 
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Figure 1. Lower bound in MeV for [rris-f mu](4GeV^) 
versus Q{GeV) from eg. ([7^. 

ting vary between 290 MeV (the up- 

per curve) and 380 Mev (the lower curve) and 
using /k = 113 MeV, Mk = 493.67 MeV. We 
take the number of active flavours rif = 3 in 
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Figure 2. Lower bound in MeV for [ms+mu] (4GeV^) 
from the quadratic inequality. 

are the lower bounds obtained by letting A^i^ 
vary |^ between 290 MeV (the upper curve) and 
380 Mev (the lower curve). Values of the quark 
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masses below the solid curves in Fig. 2 are forbid- 
den by the bounds. The horizontal lines in this 
figure are also the same lattice results as in Fig. 1. 

The quadratic bound is saturated for a i5-like 
spectral function representation of the hadronic 
continuum of states at an arbitrary position and 
with an arbitrary weight. This is certainly less 
restrictive than the extreme limit with the full 
hadronic continuum neglected, and it is there- 
fore not surprising that the quadratic bound hap- 
pens to be better than the ones from Ejv((5^) for 

= 0, 1, and 2. Notice that the quadratic bound 
in Fig. 2 is plotted at higher Q-values than the 
bound in Fig. 1. The evaluation of the corre- 
sponding QCD factor with or without inclusion 
of the 0{al) terms differs by less than 10% for 
Q2 > 4 GeV^ and we consider therefore that for 
the evaluation of the quadratic bound Q > 2 GeV 
is already a safe choice. We find that even the 
quenched lattice results of refs. |^J^ are in diffi- 
culty with these bounds. 

Similar bounds can be obtained for + 
when one considers the two-point function asso- 
ciated to the divergence of the axial current 



(17) 



The method to derive the bounds is exactly the 
same as the one discussed above and therefore 
we only show, in Fig. 3 below, the results for the 
corresponding lower bounds which we obtain from 
the quadratic inequality. We find again that the 



lattice QCD results of refs. 0-|9| for to„ -t- are 
in serious difficulties with these bounds. Notice 
that in this case the mass corrections in the QCD 
two-point function are negligible at the plotted 
Q- values. 

3. SCALAR ISOSCALAR CHANNEL 

The basic difference between the "scalar" and 
"pseudoscalar" channels is that the lowest state 
which contributes to the hadronic spectral func- 
tion in the "scalar" channel is not a pole but the 
TT — TT continuum with ,7 = and / = 0. This 
contribution provides a lower bound to the full 
spectral function: 



-lm^'(i) > 



1 



4r7i?. 



TT ■ IGtt^ V t 

■i\F(t)\^e(t-Aml), 



(18) 



where F{t) denotes the scalar-isoscalar pion form 
factor 



(19) 



with t ~ (ja + p') . In particular, the value of 
this form factor at the origin is the so called pion 
sigma term. This form factor has been well stud- 
ied in xPT at the one and two loop level (see 
e.g. refs[ pTp^ and references therein), and the 
predicted low energy shape agrees well with the 
available experimental information. With 



it is found 
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(20) 
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Figure 3. Lower bound in MeV for [md+rriu] (4GeV^) 
from the quadratic inequality. 



m 

and 



X (0.99 ± 0.02) -HO(to^) 



(r^)J = (0.59 ±0.02) fm^ 



(21) 



(22) 



The pion mass used in eq. ( pij) in ref. |22) is the 
7r+-mass. Since the electromagnetic interactions 
are neglected, we find it more appropriate to use 
the 7r°-mass instead, which in any case results in 
a lower contribution to the average quark mass. 
We also wish to point out that in the general- 
ized version of xPT (see e.g. ref. 23 1 and refer- 
ences therein,) the value of F{Q) could be sensibly 
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larger, which would result in even larger lower 
bounds for the quark masses. A value for the 
curvature of the scalar form factor is also quoted 
in ref. it involves however some extra phe- 
nomenological assumptions which go beyond xPT 
and this is why we have not included this extra 
information here. 

Using standard methods, (see ref . and ref- 
erences therein, and also refs.p5|-|27[) one can 
construct bounds on rh when something is known 
about the scalar-isoscalar pion form factor F{t). 
For this, it is convenient to map the complex q^- 
plane to the complex unit disc: 



l + z 
1 - z 



t 



Arrii- 



1 . 



(23) 



With F(0) and F'(0) as input, the resulting lower 
bound for is 



3tt Ami 



Nc 



2m^ 

Q 



1 + {4 +3-[-QV 



1 



ll£s(Qf) 

3 TT 



(24) 



Notice that this bound is somewhat similar to the 
first bound obtained from the divergence of the 
axial current in eq. (^5|) . There are however some 
interesting differences which we wish to point out. 
From the point of view of the large-iVc expan- 
sion, the r.h.s. of eq. ( p4[ ) is 1/iVc-suppressed. 
Also, from the point of view of the chiral expan- 
sion the r.h.s. of eq. ( p^ has an extra factor 
toJ [see eq. (pl[)] as compared to the r.h.s. of 
eq. (p!5|). One might expect therefore that the 
lower bound in (|2^) should be much worse than 
the one in (|l5|). Yet the bound for fh in cq.(24) 
is surprisingly competitive. The reason for this 
is twofold. On the one hand, the absence of a 
prominent narrow resonance in the J = 0, / = 
channel is compensated by a rather large contri- 
bution from the chiral loops, a phenomenological 
feature which appears to be common to all the ex- 
plored physical processes where the J = 0, / = 
channel contributes. On the other hand the fact 
that, numerically, rn^ and are not so differ- 
ent. The lower bounds for the sum of running 



masses (m„ -I- rrid) at the renormalization scale 
[ij,^ ~ A GeV^] which follow from this scalar chan- 
nel bound are the dashed curves shown in Fig 4 
below The dashed curves in Fig 4 correspond to 
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Figure 4. Lower bounds in MeV from the Scalar 
Channel for [md + m„](4 GeV^). 



the choice A^^ 

MS 
,(3) 



290 MeV (the upper dashed 
curve) and A^ = 380 MeV (the lower dashed 

'MS ^ 

curve) . The plots show that for values of Q in the 
range 1.4 GeV < Q < 1.6 GeV the lattice QCD 
results in are dangerously close to lower 

bounds obtained with very little phenomenologi- 
cal input, (just the values of the scalar form factor 
at the origin and its slope at the origin as well.) 
Notice also that the relevant two-point function 
in this channel is practically the same as in the 
pseudoscalar channel for which Q > 1.4 GeV is 
already a safe choice. 

It is still possible to improve these bounds by 
taking into account the extra information pro- 
vided by the fact that the phase of the scalar form 
factor F{t) in the elastic region 4mJ < t < IGmJ 
is precisely the 1 = 0, J~Ott — tt scatter- 
ing phase-shift and there is information on this 
phase-shift both from xPT and from experiment. 
The technology to incorporate this information 
is discussed in refs.|25-2^. We have restricted 
the phase-shift input to a conservative region: 
4mJ < t < (500 MeV)^, where a resommation 
of the chiral logarithms calculated at the two- 
loop level in xPT |^,^ is certainly expected to 
be reliable and is in good agreement with exper- 
iment. This results in improved bounds which 
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correspond to the solid curves shown in Fig. 4. 
The two curves reflect the variation within the 
quoted errors of the input parameters. Again, we 
find that the lattice determinations of refs. 
do not satisfy lower bounds which only incorpo- 
rate very little hadronic input. 

4. CONCLUSIONS 

We conclude that the lattice results of refs. 0- 
1^ are in serious difficulties with rigorous lower 
bounds which can be derived from general prop- 
erties of QCD and with a minimum of well estab- 
lished phenomenological input; and this in two 
different hadronic channels. The lower bounds 
we have derived are perfectly compatible with the 
sum rules results of refs. and earlier refer- 

ences therein. They are also compatible with a 
recent semiempirical determination of to^ made 
by the ALEPH collaboration based on the 
rate of the hadronic r-decays into strange parti- 
cles observed at LEP. 

Other work with various claims on possible val- 
ues of the ligh t q uark masses has recently ap- 
peared in refs. [p9|-^. 
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Discussions 

C. Allton, UKqcd collaboration 

Lattice theorists use a renormalization scale of 
fi — 2GeV since they don't believe perturbation 
theory down to fi ^ IGeV. 



E. de Rafael 

/ am aware of that. I understand of course that 

pQCD is only valid at sufficiently large euclidean 
scales; however, once you have pined down the 
quark masses at a given scale, you can always use 
the renormalization group equation to run them, 
down to whatever conventional scale (larger than 
^jjg) you want. 

A. Pivovcirov, Moscow 

How worse (or better?) would be your bounds if 
you switch off QCD corrections at all? 
E. de Rafael 

As I said, the pQCD corrections are large and 
positive; therefore, if ignored, the lower bounds 
for the light quark masses would be significantly 
larger. 



